Characteristics and benchmarks of entanglement of mixed states - the two qubit case 
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We propose that the entanglement of mixed states is characterised properly in terms of a proba- 
bility density function Vpi^E). There is a need for such a measure since the prevalent measures (such 
as concurrence and negativity) for two qubit systems are rough benchmarks, and not monotones of 
each other. Focussing on the two qubit states, we provide an explicit construction of 'Pp{£) and 
show that it is characterised by a set of parameters, of which concurrence is but one particular 
combination. Vpi^E) is manifestly invariant under SU{2) x SU{2) transformations. It can, in fact, 
reconstruct the state up to local operations - with the specification of at most four additional pa- 
rameters. Finally the new measure resolves the controversy regarding the role of entanglement in 
quantum computation in NMR systems. 
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I. INTRODUCTION 

Quantum entanglement is a unique resource for novel 
(nonclassical) applications such as quantum algorithms 
[l| , quantum cryptography Q , and more recently, metrol- 
ogy [3|. Thus, it plays a pivotal role in quantum informa- 
tion theory. It is also central to the study of the founda- 
tions of quantum mechanics Q • It is not surprising that 
an abiding interest in quantum entanglement persists to 
this date. 

Entanglement of a bipartite system in a pure state 
is unambiguous and well defined. In contrast, mixed 
state entanglement (MSE) is relatively poorly under- 
stood mainly because entanglement, as an observable 
(denoting a property of the state) cannot be represented 
by a linear operator in the Hilbert space. Although many 
criteria such as entanglement of formation and separabil- 
ity have been proposed, there is a realization Q that no 
single quantity can adequately represent the entangle- 
ment contained in a mixed state. It may, therefore, be 
worthwile investigating whether a complete description 
of MSE is possible in a manner such that the current 
criteria emerge as particular, albeit useful benchmarks. 

We propose in this paper a characterisation of MSE 
in terms of a suitably defined probability density for en- 
tanglement, Vp{£). The proposal is operational for any 
bipartite system. In this work, we focus on two qubit 
systems for which wc fully implement the definition. We 
find that Vp{£) is characterised by its points of non- 
analyticity of various orders which completely capture 
the information on MSE. This central result is employed 
to shed light on various aspects and manifestations of 
MSE. 

The plan of the paper is as follows. In the next sec- 
tion we review briefly the existing criteria of MSE, and 
their drawbacks. Section III proposes the new defini- 



tion in terms of a probability density function, which 
will be constructed fully for the two qubit case in section 
IV. Section V discusses several examples illustrating the 
proposal, and we show how an existing criterion - the 
concurrence- emerges as but one benchmark. In section 
VI we discuss an interesting application, viz., to the prob- 
lem of entanglement in NMR quantum computation. In 
section VII, we address the question of reconstructibil- 
ity of the state given its entanglement density. We show 
that, unlike in the case of every other definition, our pre- 
scription allows for an almost complete reconstruction of 
the state, upto local SU{2) x SU{2) operations. Section 
VIII concludes with a summary and outlook. 
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In this section we review very briefiy various criteria 
and definitions of MSE. It is not our purpose to provide 
an exhaustive description of all the definitions of MSE. 
We refer the reader to literature for details. Our in- 
tention is to merely provide a motivation and a proper 
setting for the new definition. 

Consider a bipartite spin system in a pure state, ~ 
Smi. m2 Cmim2|mim2), whcrc the expansion is under- 
stood in any separable basis. The entanglement in the 
state is unambiguously quantified by the entropy carried 
by the reduced density matrix of either of the subsys- 
tems, S[pr]. In the particular case of two qubits, equiv- 
alent criteria (in the sense of being relative monotones) 
include the degree of mixedness 1 — Trp^, the determi- 
nant \pr\, and the concurrence C = 2|c||Cj^| — cj^|C|||, 
in writing which we employ the equivalent and a conve- 
nient notation 1 (|) <-!■ ^(— ^). Operationally speaking, 
it is necessary and sufficient to measure a single quan- 
tity, the degree of polarisation Pi = P2 = P of either 
of the qubits. We note parenthetically that a knowledge 
of P allows a reconstruction of the parent state upto lo- 
cal SU{2) X SU{2) operations (LO). Equivalcntly, entan- 
glement determines the state upto LO. Indeed, writing 
IV') = cos{9/2)\ tt) +sin(6'/2)| H) in its canonical basis. 
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we obtain the relation P = |sin(?| which demonstrates 
the claim. 

A description of MSE is not that straight forward, even 
in the two qubit case. It is not difficult to realize that 
there exists no single parameter that characterises MSE 
0. Nevertheless, a number of concepts and associated 
quantities have been introduced in an attempt to cap- 
ture the entanglement content in a mixed state, the two 
prominent of them being separability 0] and entangle- 
ment of formation (EOF) 0, [1]. In addition, other con- 
cepts such as entanglement cost, distillablc entanglement, 
relative entropy of entanglement and entanglement wit- 
nesses have been introduced. It is instructive to look at 
the extent to which the above mentioned definitions (i) 
satisfy the requirements of entanglement, (iii) are quan- 
tifiable and, (iii) are equivalent. 

• Compatibility with the requirements: An entangle- 
ment measure is expected to satisfy continuity, ad- 
ditivity, subadditivity, convexity, and a nonincreas- 
ing nature under LO, together with classical com- 
munication (LOCC) Much work has been done 
in checking for the compatibility of the above mea- 
sures. It is found that (i) distillable entanglement 
violates convexity @ and that (ii) relative entropy 
of entanglement violates additivity It has not 
been established whether entanglement cost is com- 
patible with continuity, and if EOF is compatible 
with additivity In short, there seems to be 
no single measure which is consistent with all the 
above constraints. 

• Quantifiability: The measures listed above are 
quantifiable, at best, in a limited sense: Concur- 
rence which quantifies EOF is defined only for a 
2QS [ll|, and its generalisation to higher spin sys- 
tems is not available. Negativity as a measure of 
non-separability is a necessary and sufficient con- 
dition only for two qubit and qubit-qutrit systems 
0, ' higher spin systems, it is only a nec- 
essary condition. Other operational criteria such 
as majorization [T^ and reduction are, again, 
only necessary (but not sufficient) conditions for 
separability. Entanglement cost and distillable en- 
tanglement have eluded any quantification so far. 

• Mutual equivalence: All of the above criteria are 
equivalent only for a pure state. Concurrence and 
negativity are, for instance, not relative monotones 
and are hence inequivalent (l6j : States with the 
same concurrence can have differing negativities 
and vice-versa, although for any given state, its 
negativity is never greater than concurrence. 

To summarise, none of the above quantities can, by it- 
self, capture fully the entanglement that is contained in 
a mixed state. This observation suggests strongly that a 
complete description of MSE requires more than a spec- 
ification of a parameter. 



To further emphasise the need for a better descrip- 
tion of MSE, we note that a pure state description is 
almost always an idealisation and that any future exper- 
imental realization of quantum information processes will 
be with quantum systems in mixed states. An unsatis- 
factory understanding of MSE will reflect, in turn, an 
imprecise appreciation of the non-classical features that 
render quantum information processing possible. An ex- 
plicit example is provided by NMR quantum computers 
jl7l.[l8t. Here, the qubits are prepared experimentally in 
what is known as a pseudo pure state 

Pps = \{l-e)I + e\i,){i,\, 

where, 1-0) is a Bell state. The concurrence and the 
negativity of pps survive if e > i, while experimentally, 
e w 10-*^ (see 

These states would be essentially classical if we employ 
concurrence as a criterion for MSE, and no quantum gate 
operation should be possible with these states since the 
sole feature that distinguishes a quantum system from 
its counterpart is entanglement [l9[. Yet, notwithstand- 
ing the vanishing of this measure, nontrivial nonclassical 
gate operations with up to eight qubits have been re- 
ported [iOl. More recently, a 12-qubit pseudopure state 
has been reported for a weakly coupled NMR system [2l| . 
While one could entertain the possibility of QC without 
entanglement [2^ . it is perhaps more fruitful to unravel 
the sense in which the inherent entanglement - not cap- 
tured by EOF or nonseparability - is a resource for QC in 
these systems. Thus, there is a clear need to go beyond 
the above mentioned benchmarks and attempt to obtain 
a more complete description. We address this problem 
and propose an alternative definition of MSE in the next 
section. 



III. DESCRIPTION OF MSE BY A 
PROBABILITY DENSITY FUNCTION 

A. Motivation for the definition 

The new definition of MSE which we propose differs 
from the existing criteria in that we describe MSE in 
terms of a probability density for the entanglement. To 
motivate the idea, we recall that a mixed state descrip- 
tion is required when the system is an ensemble of quan- 
tum systems, each of which is in a pure state [1^. En- 
tanglement has a sharp value for each pure state, and it 
should be natural that MSE be described properly by a 
distribution defined over the microstates. 

This task is, however, not as straightforward as it 
might seem. Because of the principle of superposition, 
the ensemble description of a quantum system in a mixed 
state, as a weighted distribution over a set of pure states, 
is not unique. Expressed equivalently, there is no way of 
knowing how a system has been 'prepared', unless it is in 
a pure state for, only a pure state lip) belongs to a unique 
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one dimensional projection with an eigenvalue 1 

[131 ■ Thus, although MSE may be expected to acquire a 
statistical character, and be characterised by a suitably 
defined probability density function (PDF), care must 
be exercised such that the PDF for a given p it is not an 
artefact of its resolution in terms of any particular inco- 
herent superposition of pure states. As the first step in 
finding the way out, we consider the class of special sys- 
tems whose density operators are projection operators. 
Note that both the pure states and the fully unpolarised 
state belong to this class. 



B. Definition of PDF when p is a projection 

Consider the case p = -jjW^n where the projection 
operator Hm has rank M . Let 7i(njv/) be the subspace 
projected by 11^. Observe that for all \tl)) e TiiYij^i), 
{il)\p\il)) = const, which merely expresses the fact that 
the probability density in the M dimensional manifold 
7^(11^) is uniform. The density in the complementary 
subspace is, of course, zero. This statement is exact and 
does not depend on the eigenbasis (or any other set of 
states) chosen to expand p. The probability density for 
the entanglement associated with the subspace may be 
defined thus: 

Definition 1 (PDF when p is a projection) Let the 

state p = jj^M be a AI dimensional projection operator. 
The probability density function for entanglement of this 
state is given by 



(1) 



where dTin^^ is the volume measure for the manifold 
'Hum- 

To fix the volume measure in ([1]), we observe that the 
group of automorphisms G of the subspace 7i(n) leaves p 
invariant. The measure should naturally be invariant un- 
der this group action and is, therefore, intimately related 
to the Haar measure of G. Indeed, let H(n) be gener- 
ated by the group action on any reference state | V'o) ■ Any 
state IV') G ^(n) can be obtained by the action of some 
g G G: \ifj) = g\iJo)- Let H be the stabihzer group of 
the ray associated with the reference state. The measure 
dTiu is simply obtained by the Haar measure for G, af- 
ter factoring out the Haar measure for H [2^ . Since the 
Haar measure is invariant under the group action, and 
pure state entanglement is invariant under under LO, it 
follows that the PDF is invariant under LO. 

The extension of the pure state entanglement (one di- 
mensional projections) to states which are higher dimen- 
sional projections has thus turned out to be straightfor- 
ward and unambiguous. As we shall see in the explicit 
case of 2QMS which we study in detail, they have a rich 
structure which can nevertheless be captured by specify- 
ing a few parameters which are invariant under LO. It 



remains to further extend the definition to mixed states 
which are not projections. We take that up in the next 
subsection. 



C. PDF for any mixed state 

To extend the above definition to mixed states with- 
out any restriction, we adopt the guiding principle that 
two states which arc close to each other should possess 
"similar" entanglement densities. For example, the en- 
tanglement of a state with distinct but nearly equal eigen- 
values, should not differ from the entanglement of a com- 
pletely unpolarised system. To accomplish this, we write 
p as a weighted sum of projection operators Hm which 
satisfy the following property. Let H{IIm) be the sub- 
space (of dimA/ ) projected by Hm. We then require that 
n{UM) C H(Hm+i); il/ = l,---iV- 1, where N is the 
dimension of p. In terms of these nested projections Hm, 
we define the following: 

Definition 2 (PDF for a mixed state) Let a state p, 
be resolved in terms of nested projection operators as 
p = X]m=i '^m^m , with Hm satisfying the normalisation 
X^A/ ^A/ = 1- The probability density function (PDF) for 
the entanglement of p is given by 



N 



MS) 



ujmVum{S) 



(2) 



where the PDF for a projection is defined in ([7p. 



The definition given above is unambiguous since the 
weights can be easily determined in terms of the eigen- 
values of p. Let X\ be the eigenvalues of p, arranged in 
a nonincreasing order, belonging to the respective eigen- 
states \'ipi). The eigenstates are not unique if the eigen- 
values are degenerate, but they are of no consequence to 
us here. We first write the trivial identity 

p = (Ai - A2)Hi + (As - A3)n2 + 
• • • (Ajv-i — Xn)^n-i + AatHat 

TV 

= AmHm, (3) 

M=l 



where the projections H 



M 

satisfy the nestedness condition stated above. The 
weights ujm in ^ are easily read off as ujm = Am/Ai. 

With this identification, we see that the nonnegative 
vectors and A, defined by = (o^i, • • • , wjv) = AiA 
have natural, but rather different interpretations. The 
norm of A is a measure of the purity of the state, and lies 
in the range [1, 1/iV], the limiting cases corresponding to 
the pure and the completely mixed states respectively. 
The norm of O represents, on the other hand, the degree 
of projection onto a subspace. Thus, || Q, \\ takes its 



4 



maximum value, 1 when p is a pure projection. In any 
case, the form of p in ([3]) demonstrates the assertion made 
above, viz., that if a set of eigenvalues are close to each 
other, the state is then predominantly in the subspace 
spanned by their respective eigenstates, with only a small 
spill over to the individual states. In the other case when 
an eigenvalue is much larger than the other, the spill 
over to the projection to higher dimensional subspaces is 
small. These observations establish the physical viability 
of the definition. 

We remark that the definition of MSE is valid for any 
bipartite system, and is operational in the sense that it 
can, in principle, always be evaluated. The entanglement 
distribution is governed by the invariant Haar measure 
associated with the group of automorphisms of each sub- 
space, as also the entanglements of the pure states be- 
longing to it. Since they are invariant under LO, their 
structure cannot be arbitrary. Thus e.g., the PDF for a 
A'^ — 1 dimensional projection will be characterised by a 
single parameter - the entanglement of the pure state 
orthogonal to the subspace. Postponing an investigation 
to higher spins to a future work, we now implement the 
above definition to the most important case in quantum 
information theory, viz., the two qubit system. 



IV. PDF FOR A TWO QUBIT SPIN SYSTEM 

Two qubit systems are the most important from the 
view point of applications, and also because of the ex- 
tensive theoretical analyses that they have received. We 
focus our attention exclusively on 2QS in the rest of 
the paper. We (a) analyse entanglement in states which 
are pure projections, and (b) their extension to general 
states, (c) illustrate the distribution in a number of exam- 
ples, (d) discuss the role of concurrence, (e) the problem 
of reconstructing the state given the PDF for entangle- 
ment, and the (f) reconciliation of NMR QC with entan- 
glement. As our pure state measure, we choose concur- 
rence defined in the introduction. As pointed out, this 
choice does not amount to any loss of generality since all 
the measures of pure state entanglement are monotones 
of each other. 

We first consider the special class of states, p = 
-p-IlM; M = 1, • • - 4. The spectrum consists of only two 
eigenvalues, zero and 1/M with respective degeneracies 
4 — M and M. The two limiting cases d = 1,4 corre- 
spond to the completely polarised (pure states), and the 
completely unpolarised (mixed states) respectively. Each 
of the above cases will be analysed in detail. First the 
simplest of them all, viz., a pure state. 



A. One dimensional projections - the pure states 

The Haar measure for the case p = \4>){4>\ is trivial 
since the group of automorphisms is given by the sub- 
group consisting only of the identity clement. Thus, the 



PDF has the form 

in terms of the entanglement of The PDF has a sup- 
port only at £(/,, and the entanglement is characterised by 
a single number. Note that any other choice of pure state 
entanglement simply rescales £,j, — > £'^, in a monotonic 
manner. The form of the PDF is unaffected. It may also 
be noted that the PDF determines the one dimensional 
projection upto LO. 

B. Two dimensional projection - p — ^Tl2 

This particular class of states has the richest and the 
most interesting entanglement distribution. Since the 
definition of PDF in ^ takes care of the normalisation 
through the group volume factor, we pay no attention 
to the trace factor hence forth. The form of the 
PDF crucially depends on the nature of the subspace 
7^(112). Suppose that H{Il2) is spanned by the basis 
{\m1m2), \m1m2)}. Without any further computation, 
we see that every state tp G ^(112) is separable, giving 
a PDF which vanishes everywhere, except at f = 0. It 
is not difficult to see that the above statement holds 
for all subspaces related to the specified subspace by 
local operations, i.e., SU{2) x SU{2) transformations. 
Such an equivalence under LO is valid for other PDF as 
well. It is, there fore, necessary and sufficient to study 
PDF for 7^(112) which belong to inequivalent classes 
under LO. To that end, we construct a canonical basis 
in 7^(112) by freely employing LO. 

Canonical basis in 7^(112): Let lijj) € 7^(112). Let 
Ixi), 1X2) be orthonormal and span 7^(112). We have, 

1^) = cos y*/^\xi) + sin ^e-^^/2|x2), (4) 

where < 9 < tt, < (p < 2-k. The Haar measure is 
simply read off as dH = sin 6 dddcj). 

We assert that in any H{Il2), there is a state which is 
separable. 

Proof: The demonstration is straightforward. Let 
lxi)i 1X2) be an orthonormal basis in 7^(H2). Let the 
entanglement of |xi)j £xi ~ I sinaj. Its canonical form is 
then given by 

IXi) =cos|| TT)+sin|| ti), 

whence 

|X2)=a(-sin-|n)+cos-Ui)) (5) 

+b\ n) + c\ it) 

with the condition |ap + |6p -t- |cp 1. Let us ex- 
pand IV') in the above basis employing (jl]). The con- 
dition that £^ = yields the quadratic equation in 
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tan(|)ea;p(i(/)): 
, o sin Q 



hc)z^ — az + 



sma 



= 0. 



whose solutions are always physical, by virtue of the bi- 
jective mapping between the points on a sphere and the 
complex plane. 

Let the separable state be chosen as a basis state 
and be brought to its canonical form — (1,0,0,0)} 
in writing which we have ordered the basis states as 
{| TT), in), I it), I ii). Employing the residual LO 
(which leave |?7i) invariant), we may write the orthogonal 
basis vector as I772) = {0,x,y,z = ^1 — a;^ — y^), where 
x,y,z > 0. The subspace 7^(112) is, therefore, charac- 
terised by two non-negative parameters say, x,y. The 
PDF would also be characterised by the two parameters, 
and gets implicitly determined by ([1]). 



Determination of the PDF 
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FIG. 1: Some Typical probability density functions for 112. 
Note the solid curve, which shows all the features of 'P2{£)- 
It has a cusp at Ecusp — 0.8 and goes to zero at Smax — 0.89. 



The step function is an extreme example, where Scus 



0, 



and the other dotted curve, has Ec 



En 



The determination of the distribution function is 
rather involved, and we give the details in the appendix. 
Here, we present the results and discuss the salient fea- 
tures. 

The generic form of the PDF for 112 is shown in FIG. 
[T] (the solid curve) . We observe that it has three mark- 
ers, (i) Scusp, the entanglement at which the probabil- 
ity density diverges, as a cusp, {ii)£max, the maximum 
entanglement allowed, and (iii) V2i£max), the probabil- 
ity density at Smax- In fact, any two of them suffice to 
characterise the PDF completely. One may specify e.g., 

(Smax, 'P2i£max)), Or Cquivalcntly, {Scusp, 'P2{£raax)) for 

characterising the curve. A straightforward computation 
establishes the relations 



xy 



H = sm 



~ ^max cos fl 
1 



2.Jxy{xyEmax + z'^) 

c3/2 



(6) 
(7) 



(8) 



which allow us to determine the parameters x, y that de- 
fine 112. M is well defined by virtue of the inequality, 

'7^2 (^maa:) ^ ^l^max- 

The details of the nature of the PDF are presented in 
the Appendix, where the above statements are proved. It 
is also shown that the PDF itself is an incomplete ellip- 
tic integral. For our purposes here, it is important that 
the curve is characterised by the two locally invariant 
parameters x, y. 

More importantly, we note that unlike with the other 
measures, the state itself can be reconstructed up to LO. 
For. we can reexpress {x,y,z) in terms of the character- 
istics of the PDF thus: 



VS. 



max^ cusp 



^max COS 



cusp ) 
cusp)^ 



V(i-^— )(i + f, 



2/ — ^ (-'- ^ £max){^ ^ ^cusp) 



xy 



2 i^max 



^niax ) ( 1 
~ ^cusp^ — 



(9) 
(10) 

(11) 
(12) 



The above equation expresses the result that the entan- 
glement of a state which is a two dimensional projection 
is completely characterised by its SU{2) x SU{2) invari- 
ant parameters, which are essentially two in number. 



2. Relation with concurrence and negativity 

We briefly discuss the status of two well known bench- 
marks, concurrence and negativity in this description. 
It is, in fact, sufficient to consider concurrence since it 
bounds negativity from above. As a warm up, it is in- 
structive to look at two extreme cases which occur when 
Scusp = and Scusp = £max- In the first case, the PDF 
is a step function, terminating at some Smax- In the sec- 
ond case, the density increases monotonically, diverging 
at £max (see FIG.[T]). The relative abundance of the en- 
tangled states is more in the latter case. One may per 
se expect that the associated concurrence should also be 
larger. Interestingly, however, concurrence for a two di- 
mensional projection is related to the new parameters 
by 



C = {£„ 



^cusp) 



(13) 
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Thus, contrary to nai've expectations, concurrence - as a 
quantifier of entanglement of formation - vanishes when 
£cusp = £max- In othcr words, it is sensitive not to the 
relative abundance of the microstates at zero (or small 
entanglements) at all but, to the difference between Scusp 
and £max- In any case, C emerges as a particular bench- 
mark of the probability density, describing it only par- 
tially. 

We note that if p = lis or 114, then its concurrence 
vanishes identically. By virtue of its convexity, we con- 
clude that concurrence of any state Cp, obeys the inequal- 
ity 

Cp< (Ai-A2)Cni + (A2-A3)Cn,. 

Incidentally, entanglement distribution of a subspace 
7^(112) orthogonal to 7^(112) is the same as that of 7^(112). 
The proof of this statement is given in the Appendix. 



C. Three dimensional projection - p — ^Ils 

We now move on to the case p = whose PDF 
has a simpler structure. The simplicity is afforded by 
the fact that II3 is completely characterised by its dual, 
\ip±) -L 113. Accordingly, its PDF is characterised by a 
single parameter £±, which is the entanglement of the 
orthogonal state |'0j_). 
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FIG. 2: A Typical probability density for lis. Note the point 
of discontinuity in the derivative at £ = £± 

Let {Ixi)}? * = 2, 3, be an orthonormal basis span- 
ning the subspace Tiflla), under consideration. The 
integrating measure [26[ may be conveniently written 
as dTls = sm2P sin 26 sin^ 9 dad pd J d9, with the state 
expanded as \ip) = cos0|xi) -I- e*^""'"^) sin0cos/3|x2) — 
gi(a-7) gin gin Plxs) ■ The ranges of integration are given 
hy 9,P & [0, ^] and 01,7 G [0,7r]. Using arguments simi- 
lar to the ones employed for two dimensional subspaces, 
one may, conveniently, choose two of the basis states. 



say Xi,2j to be separable; by a suitable LO, they can be 
brought to the form | 11), \ H). 

In this basis, the state Ixs) — ci| ti) + C2I H), has the 
same entanglement, £±, 

as the state \il'±) = tJ,) — cl\ H). In this canonical 
form the state looks like: 

/i 0^ OX 

P = • I |2 

Q £iC, k|L 

V i / 

Here in terms of ci, C2, we have £j_ — 2\c\C2\. 

We have verified that the resulting probability density 
can be cast into the simple form 

Of" 1 

_^ cosh"^( — )■ (14) 

vi - £i ^> 

where = max(£,£j^). 

A typical curve for 7^3 (f) is shown in FIG. [H which 
exhibits the required characteristic. The curve possesses 
a discontinuity in its derivative at £±_ . Significantly, con- 
currence (being identically zero) fails to distinguish dif- 
ferent three dimensional projections, e.g., £_\_ = or 1, 
although their PDFs are vastly different. It simply picks 
up £ = 0, at which the probability density, in fact, van- 
ishes. 

Since 'P3{£) is characterised entirely by £±, it is clear 
that the state itself can be reconstructed upto LO. 

D. Full Hilbert Space 

Lastly, we consider the full space 7^(114), whose PDF 
is universal. This curve is obtained by using the Haar 
measure on SU{4:) [131 . Note that the curve is smooth 
everywhere, as shown in FIG.[3l 
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FIG. 3: The probability density 'P4,[£) for the entire Hilbert 
space. 
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E. A General Mixed State 

The generalisation of the pure state entanglement to 
higher dimensional projections has been accomplished so 
far. It remains to merely illustrate the nature of PDF 
for entanglement when we have a superposition of nested 
projections, as given in Q. The entanglement density, 
as defined in ([2]) docs retain information on the contri- 
bution from the constituent nested projections, with ap- 
propriate weights ujm- Indeed, each dimension produces 
a PDF with its indelible characteristic. PDF for Hi is 
highly singular, being a Dirac delta distribution. PDF 
for 112 is less singular, but has a cusp as well as a step 
function discontinuity at £max- PDF for lis is smoother, 
posessing only a discontinuous derivative at Finally, 
PDF for the fully unpolarised state, 114 is entirely smooth 
every where. Thus, the definition of mixed state entan- 
glement given in ^ captures all the features and stands 
vindicated. Recall that the weights have been so chosen 
that the continuity requirement is maintained naturally. 




0.2 0.4 0.6 0.8 



FIG. 4: The overall probability density 'P4,{£) for a typical 
mixed state, p, with eigenvalues {0.385, 0.288, 0.231, 0.096}. 
The features of the individual subspaces are vividly preserved. 
Note: The delta function is represented as a vertical line of 
height equal to its weight at the abscissa of its support 

Before going on to discuss more interesting examples, 
we pause to illustrate explicitly how 7^3 (£) may be de- 
termined, for a rather arbitrarily chosen state. Con- 
sider the density matrix p with eigenvalues {Xi} = 
{0.385,0.288,0.231,0.096} and the respective eigenvec- 
tors given by: 

IV'i) (0.998, 0.000, 0.031, 0.050) , 

IV'2) = (0.059, -0.009, -0.528, -0.847) , 

IV'3) = (0.000, 0.924, 0.325, -0.202) , 

IV'4) = (0.000, 0.383, -0.784, 0.489) . 

Now, for the 1-D subspace. we simply have £1 w 0.1, 
and the weight associated with this delta function PDF 



is (Ai — A4)/Ai « 0.25. Similarly, for the 2-D subspace 
spanned by \tpi),\il)2) it is easy to verify that the canon- 
ical basis in 112 can be chosen to be |xi) = (1,0,0,0), 
\X2) = (0,a;,y,z) = (0,0.00945,0.5290,0.8485). From 
this it follows that Emax ~ 0.9 and Ecusp ~ 0.8 using 
(111); ([7]) and ([8]). Also, we can easily see that £± w 0.6 
directly from \'ip4). 

FIG. m illustrates the PDF for this general case. Note 
how all the essential features stand out in the graph. 

V. EXAMPLES 

We proceed to study 'P{£) for states which are often 
encountered, and also those which have a natural geo- 
metric structure. The exercise serves to highlight the 
richness of the definition. 



A. Strongly separable states 

As the first example we consider states which are 
strongly separable, p = pi x p2- Conventional definitions 
of attribute no entanglement to this class of states. We 
can very readily look at what form the PDF will take by 
making simple local rotations so that pi and p2 are both 
written as ^(1 + fcicri*^^) (where i = 1,2 respectively), 
with ki,k2 > 0. Thus p acquires the form 



-diag(l + ki + k2 + kik2, 1 - fci -f fc2 - fcifc2, 
1 + ki- k2- fcifc2, 1 - ki - k2 + kik2) 

The eigenvector corresponding to the largest eigenvalue, 
in this case, is clearly | }}), therefore the 1-dimensional 
part of the PDF has no entanglement. The 2-D sub- 
space is also separable since the eigenvector with the next 
largest eigenvalue is either | }}) or | }}), both of which 
form a separable subspace with | }}). In other words, the 
PDF in 112 vanishes everywhere, except at £ = 0. There- 
fore the first nonzero contribution to the PDF comes from 
the 3-dimensional subspace. In this case too, it is the 
least possible contribution that is possible from a 3-D 
subspace because, £± = 0. Thus, our measure of entan- 
glement, the PDF, gives the minimum possible PDF to 
separable states. Thus the PDF is a simple superposi- 
tion of FIG. [Hand FIG. [31 The entanglement vanishes of 
course if the state is pure. 

B. Purely vector polarised states 

In the previous example, wc had a nonvanishing ten- 
sor polarisation which was not independent of its vec- 
tor polarisation. We now consider states which are 
purely vector polarised. These states are not factoris- 
able. Further, they are never in a pure state. Writing 
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/9 = j{l + Pi ■ &{ + p2 ■ (T2), we can bring it to the canon- 
ical form (0 = i(l + piaf -|-p20'|). An easy adaptation 
of the previous case shows again that the nonvanishing 
contribution comes from 



2. States of the form p = i(l + p ■ {o\ x (J2) ) 

We can easily utilise local rotations to align p along the 
z-axis, thus converting the density matrix to the form: 



C. Purely tensor polarised states 

These states come in three classes each of which we 
study below: 



p = -1+p[(jI ®(tI~(tI® erf) 






Vo 





1 

4 

-2ip 






2ip 
1 
4 









1. Pseudopure states 

An important, but an easily analysable state is a 
pseudo pure state which is an incoherent superposition of 
a one dimensional projection and the projection operator 
for the full space. These states are employed in NMR QC, 
and unravelling their entanglement is not without inter- 
est. Pseudopure states have the form p = + kai ■ o^). 
Unlike in the previous cases, the sign of k cannot be al- 
tered by a local transformation and it lies in the range 
— 1 < k < 1/3. The eigenvalue decomposition of p is 
given by 

1 + fc , , 1 - 3fc , 

p - ^{Htt + nu + Hb} + ——n's 



1 - e 



1 + en', 



(15) 



Thus, we have the eigenvalues {t + 2p, 



111 
41414 



with the respective eigenvectors: { ^^-^^^^-^^^ , | ||) 



I li 



lU>+iliT>.|^ From the above structure, the PDFs for 
various subspaces and the associated weights may be eas- 
ily obtained. For the one dimensional projection, we have 
the PDF and its associated weight given by 



2p 



(i + 2p) 



There is no contribution from the 2-D subspace, since 
LU2 ~ X2 ^ ^3 = 0. Considering the 3-D subspace, since 
£± = I, the probability density and the weights are read 
off as 



P3i£) = 2£; U3 



2p 



(3 + 2p) 



where e ~ — fc, and Tl]], Tl^ are the projection operators 
for the states | tT)i I ii); and are the projection 
operators for the respective Bell states 



Interestingly, uji = W3. 

Thus the states belonging to the above class simply 
have PDFs that are essentially linear with a slope varying 
from to 1, with a weighted ^-function at £ = 1. 



iv^B) = ^{in> + UT)} 



V2 



{in)-UT>}. 



The state is pure at the extremal value k = —1, and is the 
completely entangled singlet state. It is completely un- 
polarised at fc = 0; at the other extremal value k ~ 1/3, 
it is a three dimensional projection, orthogonal to the sin- 
glet state. Thus, for k < 0, £p gets a contribution from 
the Bell state ( the Dirac Delta has its support at £ = 1) 
with a weight uji — and the full space with a weight 

'^-'4 = i-sk ■ Similarly, when fc > 0, its entanglement gets 
a contribution from the full space with a weight 
and the three dimensional subspace (orthogonal to the 
singlet) with a weight The curve corresponding to 

lis is a straight line since £± ~ 1. There is no contribu- 
tion from the two dimensional projection in either case. 
We take up a discussion of the import of this example to 
NMR QC in the next section. 



States with traceless symmetric tensor polarisation 



Finally, we consider tensor polarised states in their 
most familiar - the quadrupolar - form p = ill + ' 
{a\ ® Cj)}, where Aij is a traceless symmetric matrix. 
This matrix is diagonalisable by a local SU (2) x SU (2) 
transformation. We bring the matrix to the form, A = 



diag(yl 



A 

XX 1 -^yy : 



-A. 



^yy), 



where A^x > Ayy > 0. 



In 



this basis, p acquires the form. 





A* 









A 

i + A 









where A = A^.^. 4- Ayy , and fi = A^x — Ayy . This gives us 
the eigenvalues 

(Ai,A2, A3,A4) = (i + 2A, - A + ^j, i - A-^), 
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where Ai > A2 > A3 > A4; and the corresponding eigen- 
vectors are the Bell states: 



iV^i) 




1 it)), 






lit)), 


1^3) 




1 U)), 


1^-4) 




1 ID) 



The rest of the analysis is straight forward. The proba- 
bility density function and the associated weight for the 
Bell state are simply read off as 



Vi{£) = S{£ ^ 1); Lui 



2A 



i + 2A' 



The PDF for the two dimensional projection spanned by 
\4>i 2) also has a simple expression and the weight given 
by,' 



V2i£) = 



UJ2 = 



X — fj, 
T 

4 



2A 



Since \ip4} is a Bell state, we obtain for the three dimen- 
sional subspace, 



V3{£) = 28; 



U)3 



2^ 



2A 



Finally the weight for the full space is given by ^t— TT"- 

It is remarkable that for this class of states, the details of 
the state manifest only in the weights. The density func- 
tion for each dimension is itself universal. Furthermore, 
we have constraints on the coefficients A + /i < j, and 
A>/u>0. 

FIG. 5 illustrates the PDF for a typical state of this 
class. 

Before we conclude this section, we wish to add the 
cautionary remark that although the above analysis re- 
veals entanglement in a host of states which are otherwise 
considered to be classical, it does not imply that all of 
them may be harnessed with equal facility. For instance, 
the entanglement of the uniform distribution cannot be 
accessed by standard gate operations which are unitary, 
and the corresponding PDF (FIG. [3]) has to be treated 
as a background. While more study is needed to discern 
the role of the entanglement density in various quantum 
information processes, there is one application of topical 
interest which we take up below. 



VI. ENTANGLEMENT IN NMR QUANTUM 
COMPUTATION 

As an application of the new description, we address 
the issue of the role played by entanglement in QC with 




FIG. 5; A typical quadrupolar state, with the values A = 

0. 15, jj, = 0.08. Note that the weighted delta function at £ = 

1, is not shown in the plot since the PDF goes to cx3 at this 
point. 



NMR. NMR QC employs the so called pseudopure states 



which have the form 



en' 



^ ^ I (Uni). It has been ex- 

perimentally demonstrated that the quantum logic op- 
erations used in QC are implementable with NMR, and 
we know that no quantum logic operation is possible with 
classical states. Interestingly, concurrence and negativity 
vanish when e < ^, while in experiments, e ^ 10~^. 

This has led to a debate on the role of entanglement in 
NMR QC [i^l although, as we saw, experiments clearly 
show that such states cannot be completely classical. 

The PDF constructed for the pseudopure states in the 
previous section resolves this problem naturally. First of 
all its 'Pp{£) is given by a weighted Dirac Delta which is 
nonvanishing for all e > 0, superposed on the background 
contribution from the uniform distribution. We know 
that the NMR signal is sensitive only to the pure com- 
ponent, the so called deviation deviation density matrix. 
Thus, although the uniform background is invariant un- 
der unitary operations, the one dimensional fluctuation 
is not, allowing for non-trivial gate operations. In other 
words, NMR QC exploits the excess of entangled states 
over the unpolarised background as a resource, and this 
feature is correctly captured by the PDF of the state. 
Incidentally, this analysis also raises the interesting pos- 
sibility of QC with more general pseudo projection states. 



VII. RECONSTRUCTIBILITY OF p FROM THE 

PDF 



Lastly, we return to the issue of the reconstructibility of 
the state (up to LO). We have seen that when p is a pro- 
jection, the reconstructibility is assured, by construction. 
When p is more general, the reconstruction is somewhat 
partial; we are not permitted to perform independent LO 
on various subspaces if the reconstruction is desired. In- 
deed, the action of SU{2) x SU (2) on p produces an orbit 
of dimension six, characterised by nine invariants. The 



10 



set of parameters which characterize the entanglement 
are seven in number, which may be chosen, for example, 
to be : {ni,tJ.2,tJ'3,£i,£cusp,£max,£±}- Thus we need 
two additional parameters which would determine p. To 
understand their role, we note that geometrically, 'Pp{£) 
is invariant under independent Local operations Li, act- 
ing on the subspaces 11^, where 11^ C 11;+ 1- If p is to 
be unique up to a global LO, one needs the additional 
constraint Li = U L^^^ , where L^*^ may be chosen freely. 

Let us choose Lj^"* — 1 (where 1 is the identity opera- 
tor). The nestedness condition, viz., that \ipi) g 1X2 and 

IV'4) G IIj, entails that lJ"^ and Lg*'^ get specified by two 
parameters each. In fact we need only one parameter 
each to fix the states up to discrete ambiguities. This is 
because from the 'Pi{£) we already know the entangle- 
ment of the states, which fixes one parameter. However, 
there remains a discrete ambiguity if only one of or 
is specified. 

We make the above argument more explicit. If 
we have 112 in its canonical form, it is spanned 
by Ixi) and 1x2) given respectively as: (1,0,0,0) 
and {Q,x,y,z). Therefore, we can specify = 



1 


Upper Bound on E 




Lower bound on E 


- 1 




- l-z2 




- 




/ 1 









26*0 TT 



FIG. 6: The upper and lower bounds on E from inequalities 
(|18p and (|19p are plotted above. 



Consider 112 first. Recall that we are considering 
the subspace spanned by |a) = (1,0,0,0) and = 
(0,a;,y, 2;). A general state IV') = cos |e*'^/2|a) -I- sin ||/3) 
has its entanglement given by 



E*'^/^ + 1x2) sin fe-^-^/^ by giving the values of = \2{zam{e/2)cos{e/2)e:'^/^ + xysin^{e/2))\^. (16) 



{0,(j}). Similarly, can be specified by (0_l,0j_) when 
it is expanded in the canonical basis of 112 = (1 ~ II2), 
given by |xf) = (0, 0, c/Vc^ -|- 6^, -b/^c^ + P) and 

The 



It follows from the above expression that the maximum 
entanglement allowed is given by 



\X^) = (0, Vc2 + b^,ab/Vc^ + b^, ac/Vc^ + b^)- 
above construction completes the argument. It is note- 
worthy that the question of reconstructibility cannot 
even be raised with other criteria. 



+ \J z'^ x^y'^ 



(17) 



It is further convenient to introduce the variable p, defined 
by z = Srnax^/cosp, and xy = Smax sin^(///2). It follows 
that the entanglemet of every state in the subspace scales 
linearly in £max- Therefore, we can write 



VIII. CONCLUSION 

In conclusion, we have shown that mixed state entan- 
glement has a rich structure and is properly described via 
a suitably defined probability density. We have explicitly 
implemented the definition to the most important case, 
the two qubit systems, and shown how criteria such as 
concurrence emerge as specific bench marks. Their pre- 
cise role in describing the entanglement is also clarified. 
The role of entanglement in NMR QC is resolved and, 
for the first time, the issue of reconstructibility of the 
state discussed. Nevertheless, the study is incomplete 
since possible applications to teleportation, quantum al- 
gorithms and error correcting codes are still to be ex- 
plored. The generalisation to higher spin systems would 
also provide a deeper and a better appreciation of quan- 
tum information processes. 



V{£)\ 



1 



f 



V {£/£.nax) 



il£ <£ max 

and '^('^^)l(£„„^ = if f > f max • 

We 

shall now utilise this scaling and concentrate on studying 

(l-z2)/2 



the distribution for £max 
and 



1. At £„ 



c, xy 



sin^ ( 



l-z' 



(1 - cos 61)2-1- 



1-z' 



sm ( 



(1 — cos ( 



I cos ( 



and 



V{£) 



1 

47r 



£_ 

2^ 



Jo 
Jo 



d(t)5{£^£{9,(j))) 
>S{£^-£\e,cl,)) 



IX. APPENDIX Now wc can do the integral pretty easily, and it leaves 

us with: 

The properties of the entanglement density 7^p(£) for £ sin9d0 

a two dimensional projection Ha will be worked out in '^^^> = ^ 2z (i^) sin0(l - cos0)| sin<^ol ^ 



detail here. 
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where 



COS 00 



2z (i^) sin6'(l - cos( 



However, we need this solution for cos^o to he in [—1, 1] 
Therefore: 



E < 



and £ > 



z am 4 
2 

z sin 6 — 



1- z' 



(1 — cos^ 



U{9) 



1 + 



C0S(6' + On 



(18) 



1-z^ 



(1 — COS ( 



m 



£ > 



cos{e - Bo), 9 > 29o 



(19) 



where 0q = 2 tan ^ z. 

The integral for 'P{£) now becomes: 

sin 61 ^6* 



where the integration is carried out over the region where 
inequalities (flSjl and (fT9|l are satisfied (see FIG. [6]). The 
denominator of the integrand goes to zero only at the 
boundaries, and both U{9) and L{9) have nonzero slopes 
almost everywhere (except L{9) at Oq). Therefore, if we 
look near such a point, Ot, we see that : 



{£ + L'{9,,)€ + 0{e^)f -£^ 
2£L\9b)e + 0{e^) 



thus near the points where the integrand blows up the 
behaviour is ~ l/v^i which is convergent. A special 
point to check at is when £ = 1 — z^, then near 9 = tt, 
both the terms in the denominator of the integral behave 



as ~ l/v^, which is a ~ 1/e behaviour, however, the 
sin 9 in the numerator also goes as ~ e, therefore the 
integral is convergent for this values as well. Therefore, 
we are left to consider the case when £ — z'^. In this 
case near Oq, the slope of L{9) vanishes. Therefore, for 
9 = 9o+e, £'^—L^{6) ~ e^. Thus the integrand behaves as 
~ 1/e, and we have a logarithmic divergence at £ = = 
cos /z. This is the cusp in the PDF. This will also scale 

as £max for values of £max 7^ 1- Thus £cusp = £max cos/x 

as mentioned earlier. 

It may be noted that the integral for the PDF can be 
recast into the form : 



if 5 > 



ri£) 




ii£ <z^ 



(20) 



*i ^/mt) Jt,3 V-R(t) 

where R{t) is a polynomial of degree 4 in i, if we make 
the substitution t = cos 6. This is basically an incomplete 
elliptic integral (the limits ti,t2... are obtained from the 
inequalities (HH) & (HH)). 

We lastly prove the result that the PDF for two 
complementary two dimensional projections are iden- 
tical. This follows from the fact that there is a bi- 
jective mapping from H(n2) to ^(113) which preserves 
the Haar volume (in fact this map is an SU{2) x 
SU{2) transformation). To demonstrate this we will 
once again choose our basis as |xi) and \x2) defined 
above. Now, 7^(112) consists of all states orthogonal to 
Ixi) and 1x2)- We now construct a basis for ^(nj), 
Ix'i) = (0, z/^z'^ + a;2, 0, -x/^Jz'^ + x"^) and IX2) = 



+ a;2, — V z2 + ; 



x^, zyjyj + a;2). It is easy to 
is separable and £ ' = £^^ . Furthermore, 



as 

(0, xy/\fz^ 
see that |xi 

the entanglement of a general state, ji/;) — a\xi) + P\X2)^ 
in the subspace is given by: £{a, (3) —2\— a(3z — 0^xy\. 
This is identical to! the entanglement of a\xi) + PW^)- 
Since the entanglement of each state in the subspace is 
identical to that in H(n2) and the SU{2) measure is the 
same, we will have the same PDFs in both these cases. 
In fact the SU{2) x SU{2) transformation that takes |xi) 
to Ixi), and IX2) to |x2) connects these two subspaces. 
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